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Abstract. Boxy/peanut bulges in disc galaxies have been associated to stellar 
bars. We analyse their properties in a large sample of iV-body simulations, using 
different methods to measure their strength, shape and possible asymmetry, and 
then inter-compare the results. Some of these methods can be applied to both 
simulations and observations. In particular, we seek correlations between bar 
and peanut properties, which, when applied to real galaxies, will give information 
on bars in edge-on galaxies, and on peanuts in face-on galaxies. 



1. Introduction 

Both simulations and theoretical studies have shown that bars are not verti- 
cally thin morphological features, but have a considerable vertical extent and a 
vertical structure, known as the Boxy/Peanut bulges (hereafter B/P; Combes 
& Sanders 1981, Combes et al. 1990). Comparisons between observations and 
A-body simulations have established this direct connection firmer (Athanas- 
soula 2005 and references therein). Observations have shown that both bars 
and B/P bulges are quite predominant in disc galaxies and that the correspond- 
ing frequencies are in good agreement with the link between the two structures 
(Liitticke, Dettmar & Pohlen 2000). 

We measure the peanut properties in a large sample of several hundred N- 
body simulations ran by one of us (EA) for different purposes. More information 
on these simulations and on their properties can be found in Athanassoula &c 
Misiriotis (2002) and Athanassoula (2003, 2007). In particular, we seek correla- 
tions between the properties of the bar and the properties of the B/P bulge. 

2. Methods for measuring bar and peanut strength and correlations 

In order to measure the bar strength, we use standard Fourier decomposition 
and take the amplitude of the second (m = 2) component. 
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Figure 1. Right figure: Strength of the bar (top panel) and of the B/P 
bulge (lower panel) vs. radius, for both snapshots. Left figure: Edge-on, side- 
on view of the boxy (top panel) and of the strong peanut (lower panel) bulges 
used in the right panels. 



where rrij and 6j are the mass and azimuthal angle of particle j. This can be 
done either globally or as a function of radius. In the former case the summation 
is carried out over all particles in the disc, while in the latter N s is the number 
of particles in a given cylindrical shell of radius R and C m ^ r = C m ^ r (R) . 

The vertical asymmetry and the strength of the B/P can be measured in a 
similar way, by 

N s 
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where N s can be either the number of particles in the disc component, or the 
number of particles per vertical column or cut (in which case C m ^ z = C mz (R)) 
and zq is the scale height of the initial exponential disc. In Fig. 1 we apply 
these two definitions to two simulation snapshots, one with a boxy bulge and 
the other with a peanut, or X-shaped bulge. 

We have also applied statistical methods, based on the distribution of the z 
coordinates of the particles in a given vertical cut (perpendicular to the equato- 
rial plane) of a snapshot seen edge-on, with the bar viewed side-on. We use the 
median of the absolute values and the standard deviation (g z ) to measure the 
strength of the peanut, the kurtosis to measure its shape and the skewness and 
mean to measure the strength of the buckling event i.e. the asymmetry during 
the vertical instability. 

Our methods to calculate the strength of the B/P bulge correlate nicely 
(Fig. 2, right panels) and can also be applied to real, highly inclined galaxies. 
Even more interesting, we find correlations between the strength and shape of 
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Figure 2. Left panels: Results of the statistical methods vs. radius. The 
upper three panels are for the boxy and peanut snapshots shown in Fig. 1 
and the lower one for a snapshot caught at the time of the buckling. Right 
panels: Two correlations showing how the statistical moments method agrees 
with the Fourier analysis method. 



the B/P bulge with the strength of the bar (Fig. 3). Note in this figure how some 
simulations undergo several buckling events, as previously noted by Martinez- 
Valpuesta, Shlosman & Heller (2006). 

3. Conclusions 

We presented several methods to calculate the strength of the bar and the 
strength, shape and asymmetry of the B/P bulge and found strong correla- 
tions between their results. The most important correlation relates the strength 
of the bar with the strength of the B/P bulge, the strongest bars having the 
strongest peanuts. We also find that the strength of the peanut depends on the 
number of buckling episodes it underwent, the strongest bars having undergone 
more buckling episodes (Fig. 3). Finally, we find a very interesting result about 
C2,z(R), he. about the shape of the radial density profiles along cuts perpendicu- 
lar to the equatorial plane. For strong bars, having a strong peanut or X-shaped 
bulge, this profile is more flat-topped, while for weaker bars, with more boxy- 
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Figure 3. Correlations between bar and peanut properties. Each symbol 
corresponds to one simulation. The type of symbol is related to the number of 
buckling events suffered by the bar during its evolution. Left panel: Strength 
of the B/P bulge measured with our Fourier based method vs. the strength 
of the bar. Right panel: Shape of the B/P bulge (i.e. shape of the radial 
density profiles along cuts perpendicular to the equatorial plane, measured 
by the minimum of the kurtosis) plotted as a function of bar strength. 



like bulges, it is more peaked. All the results summarized here are discussed in 
length by Athanassoula & Martinez-Valpuesta (2008, in preparation). 
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